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ABSTRACT 

The  problem  of  forecasting  Naval  shore  station  manning  require- 
ments is  examined  by  considering  the  flow  of  goods  and  services. 
The  approach  of  nrocess  analysis,  which  combines  alternative 
productive  processes,  is  used.  Linear  models  are  formulated  for 
centralized  planning  at  various  management  levels,  from  the  operation 
of  a  single  shore  station  to  that  of  a  command  composed  of  several 
such  stations.  The  application  of  the  decomposition  principle  as 
a  solution  technique  for  centralized  planning  is  further  developed 
into  a  methodology  for  decentralized  planning. 
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I.  PROBLEM  OF  MANPOWER  FORECASTING 

A.  INTRODUCTION 

An  important  problem  of  defense  management  is  to  determine  force 
levels.  The  U.  S.  budget  is  the  most  stringent  constraint  on  our 
defense  posture,  and  hence  it  is  desirable  in  the  interests  of  the 
National  security  to  efficiently  utilize  resources  at  all  levels. 

An  important  facet  of  the  Navy's  force  level  is  represented  by 
support  activities  ashore.  Since  these  activities  demand  resources 
which  would  otherwise  be  allocated  to  combat  forces,  it  is  essential 
that  such  activities  be  operated  efficiently. 

The  shore  establishment  may  be  considered  to  be  similar  to  a 
business  firm,  because  productivity  is  a  measure  of  effectiveness  of 
its  operation.  This  is  in  contrast  to  combat  forces  where  the  measure 
used  is  qualitative,  "combat  effectiveness."  Quantitive  techniques 
will  be  developed  to  determine  the  most  economical  utilization  of 
manpower  resources  within  the  bounds  of  specified  productivity  and 
personnel  restraints. 

B.  PURPOSE  AND  MODEL  REQUIREMENTS 

The  Navy's  Shore  Establishment  is  an  extremely  complex  organization 
with  many  inter-connected  components  or  subsystems.  In  order  to  be 
meaningful,  any  model  of  the  Shore  Establishment  must  represent  these 
interrelations.  It  must  impose  such  constraints  on  inputs  and  outputs 
as  are  applicable,  and  respond  to  changes  in  these  values  as  they  occur, 
Finally  it  must  be  compatible  with  the  data  that  is,  or  can  be, 
generated  to  measure  those  inputs  and  outputs. 


The  models  developed  here  will  be  concerned  only  with  those  labor 
inputs  that  can  be  considered  as  variable.  For  instance,  a  station  will 
have  one  and  only  one  Commanding  Officer  in  any  circumstance,  and  thus 
no  investigations  will  be  made  of  this  manpower  category,  and  others  which 
are  similar. 

C.  SCOPE  OF  INVESTIGATION 

This  is  an  initial  examination  of  manpower  requirement  forecasting. 
As  such  several  model  formulations  will  be  examined.  These  address  the 
problem  of  forecasting  at  different  management  levels  and  include 
centralized  and  decentralized  approaches.  Linear  formulations  are 
utilized  in  the  models  presented  here,  since  even  linear  models  may 
provide  substantial  insight  into  the  interrelations  of  complex  systems. 


II.   DEVELOPMENT  OF  THE  GENERAL  MATHEMATICAL  MODEL 

A.  STRUCTURE  OF  THE  MODEL 

The  models  formulated  for  manpower  forecasting  are  constructed  by 
using  process  analysis.  These  models  are  linear  programs  which,  when 
solved,  yield  the  least-cost  labor  mix  required  to  satisfy  the  capability 
requirements  stipulated  for  that  organization.  Process  analysis  is  an 
approach  for  analyzing  production  by  considering  alternative  technological 
processes.  Its  application  in  this  problem  is  detailed  in  Section  B 
below. 

For  these  models,  the  capability  constraints  are  as  follows: 

1.  Specified  output  levels  must  be  maintained. 

2.  Consumption  and  Production  of  intermediate  products  must  be 
related. 

3.  Upper  and  lower  bounds  on  labor  inputs  are  specified  by 
availability  of  personnel,  in  the  case  of  upper  bounds. 
Lower  bounds  more  commonly  arise  from  the  requirement 
that  the  Shore  Establishment  must  provide  sufficient 
billets  that  an  acceptable  sea-shore  rotation  can  be 
maintained. 

4.  Policy  constraints  may  arise  from  requirements  for 
specified  civilian  to  military  personnel  ratios  at  some 
stations,  or  in  certain  skill  categories. 

5.  Variables  must  be  non-negative. 

B.  PROCESS  ANALYSIS 

The  manager  of  any  organization  faces  the  problem  of  finding  the 
optimal  utilization  of  such  resources  as  are  available  in  order  to  obtain 
required  results.  Often  a  choice  must  be  made  from  among  alternative 
productive  processes  capable  of  achieving  the  desired  output.  Thus, 
the  manager  is  concerned  not  only  with  the  allocation  of  resources, 


but  also  with  the  selection  of  the  techniques  which  will  employ  these 
resources.  Process  analysis  models  these  technical  alternatives  by 
considering  that  the  choice  is  to  be  made  from  a  set  of  activities  known 
as  a  technological  matrix.  Each  activity  in  a  process  is  made  un  of  a 
combination  of  input  variables  in  fixed  quantitative  ratio  to  produce 
certain  output  variables.  Consideration  of  many  activities  permits  the 
method  to  model  the  technical  choices  available  to  the  manager.  A 
simplified  schematic  representation  of  the  application  of  process  analysis 
is  found  in  Appendix  D. 

C.  MATHEMATICAL  FORMULATION 

The  initial  formulation  is  intended  to  model  a  single  shore  activity, 
for  example,  a  Naval  Air  Station  or  Naval  Shipyard.  From  this  model 
extensions  will  be  made  to  allow  application  at  higher  command  levels. 
1 .  Notation 

Prior  to  formulation  of  the  problem  in  mathematical  terms  it 
is  necessary  to  define  the  following  notation. 

X  -  column  vector  whose  elements  (x.)  are  the  i — skill  category, 

in  units  of  manhours  per  month. 
Z  -  column  vector  whose  elements  (z.)  are  the  final  products  in 
the  i — category  (e.g.,  trained  pilots),  in  units  of  pilots 
per  month. 
Y  -  column  vector  whose  elements  (y.)  are  the  amount  of  the 
i —  intermediate  product,  (intermediate  products  are  those 
which  are  internal  to  the  station),  in  work  units  per  month. 

+■  h 

C  -  column  vector  of  costs  (c)  of  the  i —  skill  category  input 


(x.),  in  units  of  dollars  per  manhours  per  mon 


th 
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W  -  column  vector  of  activity  levels  (w.)  which  represent  the 
partial  productivities  of  the  technologies  which  make  up 
the  cost  centers. 
2.  Initial  Formulation 

Process  analysis  is  used  in  determination  of  the  flow  of  inputs, 
intermediate  products,  and  final  outputs  to  and  through  the  various  subcost 
centers. 

In  mathematical  form,  the  structure  of  the  model  set  forth  above 
becomes : 

minimize        C  X 
subject  to      Z  >_  K, 


AW  = 


K2  1  x  1  K3 

W,  X,  Y,  Z  >_  0 

where  C  and  X  are  (n  XI)  column  vectors,  (C  is  the  transpose  of  C), 
K,  is  a  column  vector  of  required  output  levels,  Kp  and  K~  are  column 
vectors  of  n  lower  and  upper  bounds  on  the  labor  category  inputs,  W 

A 

is  a  column  vector  of  m  activity  levels  at  which  the  cost  centers  will 
be  operated,  Z  is  a  column  vector  of  n  final  output  categories,  Y  is 
a  column  vector  of  n  intermediate  products,  and  A  is  a  matrix  of 
technological  coefficients  (N  X  m)  where  N  =  n  +  n  +  n  . 

The  formulation  above  has  both  X  and  W  as  unknown  vectors.  The 
formulation  of  a  linear  program  in  terms  of  the  activity  levels  of  the 
cost  centers  alone  yields  a  solution. 
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A(l) 

a^w  =  z 

A(2) 

,  such  that, 

A^   W  =  Y 

A{3) 

A^   W  =  X 

3.  Final  Formulation  as  a  Linear  Program 

In  order  to  eliminate  the  vector  of  labor  inputs  and  cast  the 
problem  in  terms  of  activity  levels  as  the  only  unknowns,  the  A  matrix 
can  be  partitioned  as  follows: 


A  = 


It  was  originally  assumed  that  intermediate  products  (Y)  will 
sum  out  to  zero  for  the  station  as  a  whole.  This  assumption  will, 
however,  lead  to  infeasibil ity ,  as  was  verified  by  empirical  computation. 
The  correct  formulation  of  intermediate  products  is  shown  in  Appendix  B 
to  be: 

A(2^  W  >_  0 

By  use  of  the  above  partition  of  the  matrix  A,  the  final  form  of 
the  linear  program  becomes  one  of  finding  the  vector  of  optimal  activity 
levels  which  will: 

minimize         CT  A^  W  -  CT  W 

subject  to  A^  W  >  K] 

A(2)  W  >  0 

K2  1  A^  w  -  K3 
W  >_  0  . 

4.  Determination  of  Optimal  Manning  Requirements 

The  above  linear  program  can  be  solved  by  any  convenient  L.  P. 

algorithm  to  yield  an  optimal  vector  of  activity  levels,  W*.  The  optimal 

manning  level  requirements  X*  can  be  determined  from: 

X*  =  a'3^  W*  . 
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III.  NATURE,  SOURCES,  AND  AVAILABILITY  OF  DATA 
FOR  MANPOWER  ALLOCATION  MODEL 

In  order  that  an  operational  model  formulation  be  useful,  it  must 
have  an  empirical  basis.  In  particular,  the  manpower  allocation  model 
delineated  above  requires  empirical  information  related  to  the  partial 
productivities  of  the  alternative  processes  which  reflect  the  various 
ways  in  which  the  cost  center  may  be  organized  and/or  operated.  Hence  it 
is  appropriate  to  examine  the  management  data  reporting  systems  currently 
in  use  in  the  Navy  to  ensure  that  the  model  is  designed  to  be  compatible 
with  available  data. 

A.  RMS  PRIME 

The  most  promising  data  source,  in  terms  of  continuing  availability, 
is  the  reporting  systems  RMS  PRIME.  This  system  has  established  a  break- 
down into  cost  centers  and  sub-cost  centers,  and  specified  work  units  to 
measure  outputs.  RMS  PRIME  reporting  will  eventually  be  extended  to 
all  stations.  The  reports  provide  the  following  data,  at  sub-cost  center 
level ,  each  month: 

1.  Civilian  and  military  manhours  and  direct  cost. 

2.  Total  manhours. 

3.  Total  costs  for  labor. 

4.  Total  work  units,  (intermediate  or  final  products) 

Due  to  the  fact  that  RMS  PRIME  is  to  be  a  Navy-wide  reporting  system, 
its  use  in  the  determination  of  required  technological  coefficients 
ensures  adaptability  of  the  model  to  any  shore  station  where  RMS  PRIME 
data  can  be  made  available. 
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B.  OTHERS 

Another  method  which  has  been  utilized  to  produce  data  required  for 
the  computation  of  technological  coefficients  is  the  industrial  engineering 
technique  of  work  sampling.  This  requires  a  work  sampling  team  which 
observes  personnel  on  the  job  and  reports  their  activities  in  various  pro- 
ductive (or  non-productive)  categories.  The  advantage  of  work  sampling 
is  that  only  actual  productive  time  is  accounted  in  that  category,  while 
in  RMS  reporting  the  tendency  is  to  report  all  time  spent  on  a  job  as 
productive  time.  The  great  expense  of  thorough  work  sampling,  however, 
outweighs  this  advantage  except  perhaps  as  a  spot-check  to  determine  the 
accuracy  of  RMS  data. 

Another  technique  which  has  been  found  to  yield  results  close  to 
those  of  work  sampling  is  that  of  questioning  cost-center  supervisors 
concerning  labor  inputs  required  to  meet  various  output  levels.  While 
less  costly  than  work  sampling  this  method  requires  interviewers  who  are 
knowledgeable  on  the  organization  and  operation  of  the  cost  centers 
involved. 


14 


IV.  APPLICATION  OF  MANPOWER  ALLOCATION  MODEL 
AT  VARIOUS  LEVELS  OF  PLANNING 

The  mathematical  model  formulated  in  Chapter  II  may  be  applied  at 
various  command  levels  within  the  overall  organization  of  the  shore 
establishment.  For  illustration,  the  model  will  be  applied  to  the  Naval 
Air  Basic  Training  Command,  which  consists  of  five  training  air  stations, 
and  nine  squadrons.  The  organization  of  this  command,  and  flow  of 
student  pilots  through  the  various  phases  of  training,  are  set  forth  in 
Appendix  C. 

The  application  of  the  allocation  model  to  a  single  training  air 
station,  termed  the  Single  Station  Model  (SSM),  is  to  be  investigated 
first. 

A.  SINGLE  STATION  MODEL 

Using  the  notation  defined  in  Chapter  II,  the  final  form  of  the 
linear  program  may  be  applied  directly  in  the  SSM.  The  linear  pro- 
gramming problem  which  must  be  solved  at  each  air  station  under  the 
command  of  CNABATRA  (Commander,  Naval  Air  Basic  Training)  is  the 
fol lowing: 

minimize  CTA^  W 

subject  to  A^W^^ 

A(2^  W  >  0 

K2  1  A^  w  1  K3 

W  >_  0 

where  K-,  is  the  vector  of  specified  levels  of  output  which  are  to  be 
maintained,  and  K~  and  K-  are  the  vectors  of  lower  and  upper  bounds  on 
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labor  inputs.  As  stated  previously,  by  solution  of  the  L.P.  for  an 

optimal  vector  of  activity  levels  W*  and  the  application  of  the  relation- 
ship, 

A^  W*  =  X* 
the  optimal  manning  level  for  the  station  in  question  can  be  determined. 

A  simplified  example  of  a  single  station  made  up  of  three  cost 
centers  is  presented  in  Appendix  D. 

B.  CENTRALIZED  PLANNING  AT  THE  CNABATRA  LEVEL 

An  extension  of  the  Single  Station  Model  above  can  be  used  to  model 
centralized  planning  (CPM)  for  the  entire  CNABATRA  organization.  Consider 
the  Naval  Air  Basic  Training  Command  to  be  made  up  of  s  air  stations. 
Due  to  the  method  of  training  pilots  in  phases,  moving  from  one  station 
to  the  next,  some  of  the  outputs  (z.)  of  the  Single  Station  Model  are 
now  inputs  to  other  stations  within  the  organization.  These  products,  e.g., 
partially  trained  pilots,  are  termed  internal  products.  In  order  to 
maintain  the  distinction,  the  vector  of  internal  products  is  designated 
T. 

Double  subscripting  is  utilized  to  distinguish  between  the  inputs, 
outputs,  etc.,  of  the  various  stations.  Thus  the  following  notation 
arises : 

X  -  column  vector  of  labor  inputs  required  by  all  stations. 

j-U 

X.-  column  vector  of  labor  inputs  at  the  j —  station,  such  that, 

J 

fX" 
Xl 

X  =  x? 


X 
s 

X. .  -  amount  of  the  i —  labor  input  required  at  the  j —  station, 
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Similarly  for  intermediate  products  (Y,  Y.,  y. .),  final  outputs 
(Z,  Z.,  z..),   cost  center  activity  levels  (W,  W.,  W.  .),  and  labor  input 
costs  (C,  C,  c..).  The  column  vector  of  internal  products  T,  is  defined 
such  that  t. .  is  the  i —  internal  product  produced  at  the  j —  station. 
It  is  necessary  to  partition  the  A^  '  matrix  of  technological  coefficients 
for  each  of  the  stations  to  reflect  the  difference  between  final  products 
and  internal  products.  Designate  the  matrix  formed  by  those  rows  of  the 
A^  '  matrix  of  the  j —  station  which  correspond  to  internal  products  as 
Aiu  such  that, 

a^'w.-t. 

J     3  3 

The  matrix  formed  by  the  remaining  rows  (those  relating  final  pro- 

(1) 


ducts  at  the  j—  station),  will  be  designated  A.  ' 


Note  that  the  A* 


matrices  of  the  Centralized  Planning  Model  differ  from  those  of  the  Single 


Station  Models  for  each  base  by  the  deletion  of  the  rows  which  make  up 

3 


aP>*. 


Finally  denote  by  A  the  matrix  of  technological  coefficients  for  the 
entire  organization.  A  is  structured  in  the  following  manner: 

Ia'1'  ! 


A  = 


a'1'* 

*(2) 
A<3' 


where 


and 


SO)" 


jflc) 


rA(k) 

rti 

o 

0 


(2)" 


A 


0 
(k) 


0  i 


0 

i(k) 


for  k  =  1,  2,  3. 
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An  additional  constraint  set  arises  from  the  fact  that  there  are  upper 
(and  possibly  lower)  bounds  on  the  total  number  of  personnel  of  each  labor 
category  which  are  available  to  the  commander  for  assignment.  The  constraints 


take  the  mathematical  form: 

s 


1  .  <  V  x.  .  <  b.   for  i  =  1 , 
i  —  £_i  U  -  i 


»  n 


or,  in  the  form  required  for  inclusion  in  the  linear  program, 


L  <  V  A(.3)  I 
-  3-1  J 


J  - 


Thus,  in  the  form  previously  standardized  the  problem  facing  the  central 
planning  commander  becomes: 


minimize 
subject  to 


CTA<3>  W 


A  W  = 


|t 

i 

Y 


or 


minimize 


Z   >_  K 

l  •    y      A(.3)  W.  < 
-J-l  J    J" 

W,  X,  Y,  Z,  T  >  0. 


cta<3>  W 

a^  W 

A<T>*  w 

a(2)  w 


L  <  y   A(3)  W.  •  B 
'  J-l       J  - 

W.j   0,  for  i=l, 


j  =  1 ,... ,  s 


18 


Note  that  the  results  of  Appendix  B  concerning  intermediate  products 
have  been  extended  to  the  formulation  of  internal  products  as  well. 

In  this  model,  knowledge  of  technological  coefficients  of  all  air 
stations  is  required  at  CNABATRA  level.  In  addition,  CNABATRA,  upon 
finding  the  optimal  solution  to  the  linear  program,  assigns  personnel  to 
the  stations  in  his  command  and  specifies  the  cost  center  in  which  they 
will  be  placed,  as  well  as  the  activity  level  at  which  that  cost  center 
will  be  operated. 

A  simplified  example  of  Centralized  Planning  applied  to  CNABATRA  is 
presented  in  Appendix  E. 
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V.  COMPUTATIONAL  ASPECTS 

A.  SOLUTION  OF  THE  LINEAR  PROGRAM 

The  extension  of  the  Single  Station  Model  formulation  to  centralized 
planning  at  the  CNABATRA  level  produces  a  linear  program  which,  though 
large,  is  within  the  capabilities  of  computer  software  packages  such  as 
MPS/360  for  solution.  For  example,  within  the  CNABATRA  organization  an 
average  station  could  have  approximately  200  labor  inputs,  50  inter- 
mediate and  final  products,  and  30  to  60  cost  centers. 

With  the  capability  to  handle  4095  constraint  equations  and  unlimited 
variables,  MPS/360  can  be  programmed  for  the  centralized  planning  model 
of  CNABATRA  and  its  eight  component  stations.  Note,  however,  that  at 
command  levels  above  that  of  CNABATRA  a  centralized  planning  model  quickly 
becomes  too  large  to  be  handled  by  any  existing  computer  program. 

B.  DECOMPOSITION  PRINCIPLE  AS  A  SOLUTION  TECHNIOUE 

The  principle  of  decomposition  [DANZIG,  1963]  was  introduced  as  a 
solution  technique  for  large  linear  programs  representing  complex  super- 
systems  in  which  the  sub-systems  operate  almost  independently  but  are 
tied  together  by  some  set  of  constraints  on  the  super-system. 

The  general  form  of  the  constraint  set  of  a  problem  readily  amenable 
to  the  decomposition  principle  is  as  follows: 

Ai  xi  ■  bi 

A2  X2  =  b2 

mm        m 

A*  X,+  A*  X0.   .   .+A*  X      =  b  .. 
112  2        mm        m+1 

X.   0  for  all  j 
,1  - 
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* 
Where  the  A.  and  the  A.  are  matrices,  the  X.  are  vectors  of  activity 

levels,  and  the  b.  are  vectors  of  constants. 

Inspection  of  the  A  matrix  of  the  Centralized  Planning  Model  indicates 

that  the  CPM  may  indeed  be  put  in  this  form,  that  is: 

minimize  7  C.  a\   W. 

subject  to  A,  W, 

A2  W2 


AW  =  B 

s  s 


2     2  s 


* 

W  =  T 


L  <  A,(3)  W,  +  A^3)  W9  +  ...  +  A^3)  Wc 
—  2     2         s     J 


1 

(1) 

2 

(2) 

s 

(3) 

(4) 

(5) 

W.  >_  0  for  j=l ,. . .  ,s 


where 

B. 
3 


'I.1 

j  J 

iV 
ft  J 


Constraint  sets  (1),  (2)  and  (3)  are  independent  of  each  other,  that 
is,  the  stations  in  CNABATRA  are  only  jointly  constrained  by  (4)  and  (5), 


21 


VI.  DECENTRALIZED  PLANNING  AT  CNABATRA  LEVEL 

The  basic  philosophy  of  decentralized  planning  is  to  enable  a  super- 
system  planner  to  optimize  the  operation  of  the  entire  organization  under 
his  cognizance  without  detailed  information  and  central  direction  for  each 
decision.  The  planner,  who  must  satisfy  overall  constraints  which  bind 
together  the  entire  organization,  simply  specifies  the  prices  which  sub- 
systems must  pay  for  scarce  resources.  By  adjustment  of  those  prices  he 
"forces"  the  subsystems  to  seek  resource  requirement  solutions  which 
yield  an  optimal  solution  to  the  supersystem  problem.  Subsystem  managers, 
who  have  the  detailed  knowledge  of  the  productive  processes  necessary  to 
determine  optimal  technological  alternatives,  are  required  to  adjust  their 
operations  in  the  face  of  resource  price  changes  in  order  to  optimize  the 
cost  of  operation  at  the  new  prices. 

In  this  way,  the  supersystem  planner  is  concerned  only  with  the 
generation  of  resource  requirements  which  are  a  feasible  solution  to 
supersystem  constraints.  Decisions  concerning  technological  process  and 
skill  category  mixers  are  delegated  to  the  subsystem  manager  closest  to 
those  details. 

The  above  considerations  may  be  applied  to  obtain  a  decentralized 
planning  model  for  the  CNABATRA  organization.  That  is,  one  in  which  the 
CNABATRA  Planning  Staff  can  be  assured  of  reaching  an  optimal  least-cost- 
labor  solution  to  the  Training  Command,  without  being  concerned  with  the 
organization  or  operation  of  any  of  the  air  stations  which  make  up  the 
Training  Command. 

It  was  shown  in  Chapter  V  that  the  only  elements  of  the  entire  constraint 
set  which  jointly  affect  all  the  stations  are: 
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and 


A^*  W,  +  Ap>*  W2  +   ...+A<1)*Ws       .         T  (1) 


LiAf'vf1  »2t...Hf»3    i        B  t« 


Constraint  set  (1)  above  relates  internal  products.  These  internal 
products  are  final  products  in  the  Single  Station  Model  and  were  redes- 
ignated as  "internal"  because  they  form  inputs  to  stations  at  some 
advanced  point  in  the  training  syllabus. 

In  the  case  of  the  Naval  Air  Basic  Training  Command  these  internal 
products  are  partially-trained  student  pilots  advancing  from  one  phase  of 
training  to  the  next.  Thus  the  final  output  of  completely  trained  pilots 
is  the  number  of  students  who  entered  the  program  less  the  attrition 
at  each  phase  of  training.  Sufficient  data  exist  to  determine  the  mean 
required  output  at  each  phase  of  the  training  syllabus  in  order  that  the 
final  output  of  the  system  meet  the  specified  level.  In  this  manner 
output  levels  for  each  training  station  may  be  specified  as  elements  of 

en*  (i) 

K,,  and  the  Av  '     constraint  set  returned  to  Av  .  In  this  way,  deter- 
mination of  the  necessary  levels  of  internal  products  allows  them  to  be 
considered  as  final  products  again  and  constrained  in  the  same  manner  as 
in  the  original  Single  Station  Model.  The  station  which  received  internal 
product  (t.  .)  will  now  consider  that  quantity  as  a  labor  input  which  has 
as  a  lower  bound  the  level  specified  as  an  output  requirement  to  the 
producing  station.  In  this  way  constraint  set  (1)  is  eliminated  as  one 
held  in  common  by  all  stations. 

Constraint  set  (2)  relates  overall  personnel  availabilities  and  is 
as  follows: 

L  <  A^  W,  +  a!,3)  W0  +  ...  +  A^3)  W  <  B 
—  c         c  s    s  — 


23 


From  Chapter  II ; 


A^  W.  =  X. 
J    J     J 

for  j  =  1 , . . . ,  s 


which  yields 


L  <  X,  +  X0  +  ...  +  X  <  B 
—      c  s  — 


Note  that  in  comparison  with  the  general  form  of  the  decomposition 
algorithm  the  A.  are  identity  matrices  in  this  formulation.  (See  Appendix 
F,  equation  (3)  ) . 

A  supersystem  linear  program  may  now  be  written: 

minimize 


CT  X 


subject  to 


L  <  X-.  +  X0  +  ...  +  X  < 
-      2        s  — 


(3) 
(4) 


and  further  constrained  that  the  X.  satisfy 

3  J 


A]  W1 


Ap  W2 


=  B. 


(5) 
(6) 


AW    =  B 
s  s     s 

W.  i  0;  for  j  =  1 ,  ...,s. 


(7) 


where  again, 


Z. 
J 


Thus,  the  decentralized  Planner's  concern  is  that  the  vector  sum  of 


the  solutions  X.,  submitted  by  the  stations  in  his  command,  represent  an 

* 
optimal  feasible  solution  to  (3)  and  (4)  above.  X.  is  determined  by 

th  * 

the  j —  station  by  finding  the  optimal  solution  W .  to  the  subprogram: 


minimize 


cTa<3»w. 

J  J    J 
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subject  to  AVMrW.  >  K,. 

a(2)  w.  >  0 
J       J  - 

K0.    <   A^   W.    ■:   K0, 

W.  >  0 
J  - 

and  solving  for  X.,  by  use  of  l\\   '   W.  =  X. 

1 .  Decomposition  Principle  Applied  to  Decentralized  Planning 

As  set  forth  in  Appendix  F,  the  Decentralized  Planner  "forces" 
the  station  to  submit  manpower  requests  which  yield  a  feasible  solution  to 
the  overall  problem,  by  changing  the  prices  to  be  charged  for  labor  inputs. 
That  is,  in  response  to  an  increase  in  the  price  of  labor  inputs  which 
violate  the  upper  bound,  stations  will  seek  solutions  which  require  less 
of  that  input.  Similarly,  a  decrease  in  the  price  of  a  labor  input  will 
cause  stations  to  attempt  to  utilize  more  of  that  input,  with  the  result 
that  the  lower  bound  may  be  satisfied.  It  may  be  necessary  that  several 
iterations  are  required  before  feasibility  of  the  supersystem  program 
is  reached. 

2.  Initiating  the  Algorithm 

Assume  that  requests  for  manpower  allocation  X?,  x2,...,  X°, 
have  been  submitted  by  the  stations  as  optimal  solutions  to  their  respective 
subprograms  using  the  actual  wages  paid  in  the  different  labor  categories. 
Further  assume  that  these  X°  do  not  form  a  feasible  solution  to  the  Master 
Program,  i.e.,  some  upper  or  lower  bounds  are  violated.  The  phase  I-phase 
II  method  of  the  Revised  Simplex  Method  will,  in  general,  be  required  to 
initiate  the  algorithm. 
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CNABATRA,  as  a  Decentralized  Planner,  must  change  prices  on  labor  inputs 
in  such  a  way  that,  given  the  new  prices,  the  stations  will  change  their 
optimal  labor  mix  to  yield  an  optimal  solution  to  the  Master  Program. 

To  determine  the  new  set  of  prices,  CNABATRA  must  solve  this  L.P., 
termed  the  restricted  master  program. 

minimize 

I       CT  X°  A . 
j  =  l   J   J   J 


subject  to 


L  <_  X°  x]  +  X°  x2  +  ...  +  X°  As  <  B  (n  constraints) 

(8) 


x1  =1 

Ap  =  1       s  constraints 


(9) 


=  1 


Let  the  lower  bound  L  be  equal  to  zero,  for  this  example.  The 
constraints  on  the  subprogram  ensure  that  X.  =  0.  Note  that  the  values 
of  the  a.  are  all  one  because  there  is  only  one  basic  feasible  solution 
to  the  subprograms  which  is  know. 

Assuming  that  the  vector  X.  has  n  elements,  the  dual  of  this  Master 

Program  will  have  n  +  s  dual  variables.  The  value  of  each  of  these 

x 

dual  variables  can  be  thought  of  as  the  opportunity  cost  (or  shadow  price) 

2 
of  the  input  constraint  with  which  it  is  associated.   Let  the  row  vector 

of  dual  variables  be  designated: 

V  =   [V],  v2,...,  vn  +s  ] 


2 
Hadley,  G.,  Linear  Programming,  Addison-Wesley ,  p.  484,  1963, 
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In  the  standard  Simplex  minimization  algorithm,  a  vector  solution  not 
presently  in  the  basis  will  be  considered  for  entry  into  the  basis  if  its 
"Z.  -  C"  is  positive,  i.e.,  its  indirect  cost  is  greater  than  its  direct 
cost.  That  vector  X.  will  be  selected  for  entry  which  corresponds  to  the 

max  (Z.  -  C).  In  the  restricted  master  program  only  those  vectors  which 

j    J    J 
form  the  basis  appear.  Thus,  the  "Z.  -  C."  for  vectors  not  in  the  basis 

cannot  be  calculated  because  those  vectors  have  not  even  been  determined. 
The  decomposition  principle  provides  a  method  for  generation  of  new  non- 
basic  vectors  only  if  they  are  likely  candidates  for  entry  into  the  basis, 

The  indirect  cost,  Z.,  is  calculated  in  the  Simplex  method  by  the 
following  equation: 


wh 


ere  Cp  is  the  transpose  of  the  column  vector  of  costs  of  variables  in 


,-1  • 


the  basis,  B   is  the  inverse  of  the  present  basis,  and  a.  is  the  column 
vector  of  coefficients  of  a.,  j  not  in  the  basis. 


From  Duality: 


V  =  C 


T  n-l 


thus 


Z.  =  V  a. 
3  3 


The  vector  a.  of  the  Standard  Simplex  method  corresponds  in  this  case 


to;  (refer  to  constraints  (8)  and  (9)) 

,o 


a . 
3 


where  e^  is  a  s-dimensional 
column  vector  with  a  1   in  the 
j —  position  and  zeros  else- 
where. 


Then  the  Z.  for  the  j—    station  is 


Z.     =     v,x-,  .  +  v0x0.  +   ...  +  v       x      .  +  v 
3  1    lj         2  2j  nx     V         nx  J 
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or: 

Z.  -  C.  =  V  X°  -  CT  X°  +  v  .  .  =  (V  -  CT)  X°  +  v  A. 
J    J       J    J  J    nx+j         y     j    nx+j 

The  costs  and  the  values  of  the  dual  variables  are  all  known.  A  new 
vector  to  be  brought  into  the  basis  must  satisfy: 

max  [(V  -  CT)  X  +  v  ..]  >  0 
3  nx+jJ 

Note  that  v  ,.  is  a  constant,  thus  the  max  I'  (V  -  C.)  X.]  must 
n  +J  J   J 

th 
be  determined  by  the  j —  subprogram  to  decide  the  vector  which  is  to  be 

entered  into  the  basis.  It  is  desirable  to  maintain  the  form  of  the  sub- 
program as  a  minimization  problem,  in  order  to  preserve  the  illusion 
(at  the  station  level)  that  "prices"  are  simply  being  changed.  This  is 
accomplished  by  changing  the  objective  function  to  the  subprograms  to: 
min  [(cj  -  V)  Xj] 
In  terms  of  activity  levels,  the  subprogram  problem  for  the  j — 


station  becomes 


minimize      (CT  -  V)  A^  W 


subject  to    Av'  W. 


A(.2)  W.  >  0 


J    J 
J  -  Ij 

.  > 

J  - 


K0.  <  A^  W.  <  K0. 


W.  >  0 
J  - 


After  solution  of  the  subprograms  with  the  new  prices,  CNABATRA  will 

receive  new  manning  level  requests,  X,  ,  X?,  ...  ,  X  ,  and  corresponding 

costs  of  operating  each  station,  these  artificial  "operating  costs" 

are  termed  C .. 
J 
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The  negative  of  the  C.  are  added  to  v  +.s  and  the  X.  corresponding 

J  Y 

to  the  max  [-C.  +  v  +.]  >  0  is  selected  for  entry  into  the  basis. 
Consider  for  example  that: 

maX  ["CJ  +  Vn/j]  =  ["Ct  +  \^ 

Then  xi  will  enter  the  basis  of  the  restricted  master  program  as 
the  coefficient  of  a  new  variable  x..  Thus  the  restricted  master  has 
the  form: 

minimize 

clx?x,  +  clx2x«  +  ...  +  cl  (X°A.  +  xlxl)  +  ...  +  CTX°A 

CCC  ttt  t    t  SSS 


subject  to 


X?A,  +  X°a0  +  ...  +  X°a.  +  xJa!  +  ...  +  X°a     <  B 

C    C  tt  tt  SS   — 


x2  =  1 


(10) 


v  4  =  ] 


=  1 

Note  that  at  each  iteration  only  one  new  vector  of  inputs  and  one 
corresponding  variable,  a,  is  introduced  into  the  restricted  master. 

3.  Subsequent  Iterations 

Solution  of  the  new  restricted  master  leads  to  another  vector  of 
dual  variables  and  computation  of  a  new  set  of  "prices".  The  prices  are 
specified  to  the  stations  and  a  new  determination  of  minimum  operating  is 
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made,  allowing  entry  of  a  new  vector  solution  into  the  basis  of  the  master 

This  iterative  process  is  continued  until  there  is  no  j  such  that 

-C  +  v  ..  >  0.  This  condition  ensures  optimal ity  in  the  master  program. 

In  the  general  case  where  the  lower  bounds  are  greater  than  zero, 
the  only  modification  to  the  above  is  that  the  constraint  set  (8)  becomes 
twice  as  large. 


That  is 


L  <  X?x  ,+X?  x0+   ...  +  X°  A  < 
—  1     C      c  s  s  — 


becomes 


X°  yx°  y  ...  U»»scB 

A  vx2  v  •••  +  Xs°  h-1 

Thus  there  2n  +  s  dual  variables,  the  algorithm  is  operated  as 

A 

indicated  in  any  case. 

B.   IMPLEMENTATION 

The  decentralized  planning  model  could  be  implemented  in  the  following 
manner  at  CNABATRA. 

CNABATRA,  having  specified  output  levels  which  must  be  maintained 
by  each  of  the  stations  in  his  command,  directs  that  each  station 
commander  submit  proposed  manning  levels  which  will  produce  the  required 
outputs  at  a  minimum  expenditure  for  labor  inputs,  considering  actual  wage 
levels.  At  CNABATRA  the  proposals  are  combined  and  a  restricted  master 
program  constructed  using  the  proposed  levels  as  the  X.  vectors  of  the 
constraints.  CNABATRA  staff  performs  the  computations  to  determine  the 
necessary  values  which  are  specified  to  the  stations  as  "new  prices"  of 
the  labor  inputs,  and  the  stations  are  required  to  propose  new  manning 
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levels  at  the  new  prices.  A  new  vector  is  selected  for  entry  into  the 
restricted  master  basis,  a  second  set  of  artificial  "prices"  determined 
and  new  solutions  requested  from  the  station  subprograms.  This  iterative 
procedure  is  continued  until  optimal ity  is  reached  in  the  master  program. 
The  flow  of  information  can  be  represented  as  follows: 


-4 


CNABATRA 
PLANNER 


Terminates  when 
optimal ity  reached 


Submits  new 
i  manning  level 
proposals 


Specifies  new 
prices  for  labor 
inputs  


TRAINING 
AIR  STATIONS 
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VII.  OTHER  PLANNING  INFORMATION  AVAILABLE 
FROM  LINEAR  PROGRAM  SOLUTION 

The  utilization  of  the  IBM  linear  programming  package,  MPS-360, 
greatly  enhances  the  value  of  the  models.  Its  particular  advantage 
arises  in  the  use  of  the  Decentralized  Planning  Model.  Upon  solution  of 
the  restricted  master,  values  of  the  dual  variable  activity  can  be  putputted 
and  a  quick  multiplication  will  yield  the  new  objective  function  to  be 
specified  to  the  subprogram. 

An  advantage  of  MPS-360  is  the  ease  with  which  post-optimal ity 
analysis  of  the  solution  thus  generated  can  be  performed. 

The  following  investigations  can  be  done  by  the  system: 

1.  Over  what  range  of  values  can  the  costs  of  labor  inputs 
be  varied  without  affecting  the  optimal ity  of  the  present 
solution? 

2.  Over  what  range  of  values  may  the  "right  hand  sides"  be 
varied  before  feasibility  of  the  present  solution  is 
affected?  For  instance,  how  low  may  the  upper  bound  on 
personnel  of  a  given  labor  category  be  set  before  a 
different  solution  is  required? 

3.  Over  what  range  may  the  values  of  the  technological 
coefficients  be  varied  without  affecting  the  optimal ity 
of  the  present  solution?  This  may  be  used  to  determine 
the  impact  of  an  anticipated  change  in  technology  of 
the  present  operations  of  an  air  station  or  a  single 
cost  center. 

4.  Use  can  be  made  of  parametric  programming  to  answer 
questions  which  are  an  extension  of  those  posed  above. 
For  instance,  if  it  is  known  that  the  required  Pilot 
Training  Rate  can  be  expected  to  vary  over  a  given 
range  in  subsequent  time  periods,  parametric  program- 
ming can  be  utilized  to  determine  if  the  present 
solution  will  remain  optimal  over  the  entire  ranqe, 

or  some  portion  of  it,  and  can  suggest  a  solution  that 
will  come  "closest"  in  some  sense,  to  an  optimal 
long-run  manning  posture. 
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Both  the  post-optimal ity  analysis  and  parametric 
programming  techniques  can  be  applied  to  as  many 
individual  values  as  desired,  as  well  as  to  entire 
rows  and/or  columns,  allowing  the  user  an  almost 
unlimited  number  of  situations  which  can  be  simulated, 
and,  in  many  cases,  problems  that  can  be  foreseen  and 
avoided. 
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VIII.  CONCLUSIONS 

The  quantitative  aspects  of  forecasting  manpower  requirements  for  a 
Naval  Shore  Activity  have  been  examined  and  mathematical  models  have  been 
formulated.  Based  on  this  study,  the  following  conclusions  have  been 
reached: 

1.  The  activities  of  the  Naval  shore  establishment  may  be 
modelled  using  process  analysis  and  linear  programs 
developed  to  yield  least-cost  solutions  for  manning 
levels  required  to  meet  given  output  specifications. 

2.  In  the  centralized  planning  model,  detailed  information 
on  each  component  activity  must  be  available  to  the 
central  planning  staff.  In  contrast,  for  decentralized 
planning  each  subsystem  manager  is  responsible  for  the 
efficient  operation  of  activities  under  his  control. 
The  overall  planner  coordinates  these  activities  by 
consideration  of  constraints  known  only  at  his  level. 

3.  Decentralized  planning  by  decomposition  may  be  shown 
to  reach  the  same  optimal  allocation  solution  as 
centralized  planning. 

4.  The  decentralized  planning  technique  requires  that 
the  cost-center  manager  (who  is  closest  to  the  actual 
operations)  make  organizational  decisions  utilizing 
his  experience  and  encouraging  his  ingenuity. 

5.  Decentralized  planning  can  be  utilized  at  any  manage- 
ment level,  even  at  such  low  levels  that  no  model  is 
required  for  organization  below. 
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APPENDIX  A 

General  Assumptions  of  Process  Analysis  Model 

A.  Explicit  Economic  Assumptions  of  the  Models  are  as  follows: 

1.  Axiom  of  proportionality  -  requires  that  the  partial  productivities 
of  all  inputs  are  independent  of  activity  levels.  That  is,  the  tech- 
nological coefficients  of  the  cost  centers  are  not  a  function  of  the 
activity  level  at  which  the  cost  center  is  operated,  or  "constant  returns 
to  scale." 

2.  Axiom  of  additivity  -  requires  that  the  combined  output  of  two 
technological  processes  operated  simultaneously  is  equal  to  the  sum  of 
the  individual  outputs  when  operated  at  the  same  activity  levels  that 
is,  technological  coefficients  are  not  functions  of  interaction  with 
other  cost  centers. 

3.  Maxim  of  economic  efficiency  -  requires  that  whatever  activity 
levels  are  set  cost  centers  will  each  be  operated  as  efficiently  as 
possible. 

4.  Capital  investment  held  constant  during  the  period  for  which 
requirements  are  to  be  predicted.  In  the  long  run  such  changes  in 
capital  investment  will,  however  be  reflected  in  changes  in  technological 
coefficients  to  reflect  increased  (decreased)  efficiency  of  operation. 
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APPENDIX  B 

FORMULATION  OF  INTERMEDIATE  PRODUCT  CONSTRAINTS 

"A  mathematical  problem  which  is  to  correspond  to  physical  reality 
should  satisfy  the  following  basic  requirements: 

1.  The  solution  must  exist. 

2.  The  solution  should  be  uniquely  determined. 

3.  The  solution  should  depend  continuously  on  the  data, 
(requirement  of  stability). 

Any  problem  which  satisfies  our  three  requirements  will  be  called  a 
properly-posed  problem." 

In  order  that  the  formulation  of  the  allocation  model  as  a  linear 
program  be  "properly-posed,"  the  requirements  stated  above  must  be 
restated  as  follows: 

1.  A  feasible  solution  must  exist. 

2.  While  alternative  optimal  labor  mixes  may  exist,  the  cost 
of  the  optimal  solution  will  be  unique. 

3.  At  the  least,  piece-wise  linearity  of  outputs  as  a  function 
of  inputs  should  be  realized. 

The  assumption  that  intermediate  products  sum  to  zero  has  intuitive 

appeal,  as  the  intermediate  products  generated  as  a  positive  output  from 

one  cost  center  will  be  inputs  (negative)  to  another  cost  center.  It 

will  be  shown  that  this  formulation  results  in  an  infeasible  problem. 


Courant,  R.  and  Hilbert  D.,  Methods  of  Mathematical  Physics,  Vol.  II, 
p.  227,  Interscience,  1962. 
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A.  DISCUSSION 

The  L.P.  model  has  been  formulated  as: 

minimize  C  A   W 

subject  to  A^  W  >  K] 

A(2)  W  =  Y 
K2  1  A^  w  1  K3 
W  >_  0 

The  usual  assumption  made  is  that  the  intermediate  products  sum  to 
zero  throughout  the  station,  which  would  yield 

A^  W  =  0 
Consider  first  the  case  of  one  technology  per  cost  center.  Thus  if 
there  are  N  cost  centers,  W  is  a  N  x  1  column  vector. 
Let:  f  =  number  of  final  products 

I  =  number  of  Intermediate  Products 
N  =  number  of  Cost  Centers 

Theorem  1 :  I  •  N  -  f 

Proof:  Each  cost  center  has  some  type  of  unique  quantifiable  output, 
otherwise  the  elements  of  the  A  matrix  for  that  cost  center  would  be 
undetermined.  The  least  number  of  intermediate  products  correspond  to 
the  case  where  each  cost  center  produces  one,  and  only  one,  output, 
(either  final  or  intermediate).  Thus  the  minimum  I  is  equal  to  N  -  f. 

In  general,  I  will  be  greater  than  N  -  f. 

(2)  (2) 

A         is  an  I  x  N  matrix,  since  there  is  one  row  of  Av  '  corresponding 

to  each  intermediate  product,  and  one  column  corresponding  to  each 

cost  center 


37 


(2) 
Thus  the  maximum  rank  of  the  Av  '   matrix  is: 

r[A^]  =  min  [I,  N] 
(1) 


(1) 


Theorem  2:  Given  that  the  Av  '   matrix  is  of  full  rank,  f;  The  rank  of 

f  2) 
Av  '  must  be  less  than  or  equal  to  (N  -  f)  for  a  feasible  solution  to 

exist. 


f.  Consider  the  following  two  cases 
(2) 


Proof:  By  Theorem  1 ,  I 

Case  1:  I  =  N  -  f 

By  (1)  above,  the  maximum  rank  of  hKC}   is  the  min  [I,  N],  and  the 
theorem  is  satisfied. 

Case  2:  I  >  N  -  f 

By  contradiction,  assume  that  r(l\KCI)    •  N  -  f  and  a  feasible 
solution  exists.  Consider  the  matrix: 

(1) 


(2) 


A  = 


A 


(2) 
(2) 


The  assumption  above  that  the  r(Av  ;)  =  M  >  N  -  f,  implies  f  linearly 
independent  rows  of  A^  '   and  M  linearly  independent  rows  of  A^  '. 
Denote  the  f  rows  from  A^  '  as  j  B-  H       and  the  M  linearly  independent 
rows  of  A^  as  )   B?2^  . 


th 


The  members  of  the  set  i$:      >      cannot  be  linearly  dependent  on 
ose  of  )   B-  ' { 


Otherwise  that  member  could  be  written: 


(1) 


3-1  J  J 


(2) 


And  multiplying  on  the  right  by  W 


41'  w 


=  l      X  .  B(.2)  W 

3=1  J   J 


XB(2»  H 


=  A-  0  =  0 


where  x  =  [^  ^,  ...  ,x  m] 
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That  is,  3-   W  =  0,  which  does  not  satisfy  the  constraints  on  final 
outputs.  Thus  there  can  be  no  linear  dependencies  between  j  g:  c 
and  j   3.  \    .     Hence  the  rank  of  A  is  equal  toM+f>N-f+f=N, 

and  the  contradiction  is  reached  r  (A)  >  N. 

(2) 
Thus,  rank  of  l\     '   must  be  less  than  or  equal  to  n  -  f. 

B.  CONCLUSIONS 

By  Theorem  2,  no  feasible  solution  exists  to  the  linear  program  with 

(2)  (2) 

the  constraint  set  Av  '   W  =  0  unless  the  rank  of  Pc   '   is  less  than  or 

equal  to  N  -  f.  Since  it  is  unlikely  that  actual  linear  dependencies 

(2) 
between  rows  of  the  Av  '  matrix  would  be  detected  due  to  the  imprecise 

nature  of  the  derivation  of  technological  coefficients  and  computer 

round-off.  In  consideration  of  the  above,  the  proper  formulation  of 

the  linear  program  is  as  follows: 

T  (3) 
minimize  C  kK    '   W 

subject  to  A^  W  >  K] 

A^2)  W  >  0 
K2  1  A^  w  1  K3 
W  >  0. 
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APPENDIX  C 
SYNOPSIS  OF  CNABATRA  ORGANIZATION 

The  Naval  Air  Training  Command  consists  of  five  airfields,  with  a 
total  of  nine  training  squadrons.  The  flow  of  pilots  through  the  command 
is  as  follows.  Primary  flight  training  is  conducted  by  VT  -  1  at  Saufley 
Field  for  all  student  pilots.  Upon  completion  of  the  primary  syllabus 
those  students  selected  for  jet  training  advance  to  Basic  Jet  Training  at 
Meridian,  Mississippi,  (VT  -  7  and  VT  -  9) ,  then  to  Sherman  Field  (VT  -  4) 
for  aircraft  carrier  qualification  and  finally  are  sent  to  either  Kingsville 
or  Chase,  Texas  for  advanced  jet  training. 

Students  selected  for  propellor  aircraft  training  upon  completion  of 
the  basic  syllabus  advance  to  Whiting  Field,  Milton,  Florida  for  basic 
prop  training  (VT  -  2  and  VT  -  3),  then  return  to  Saufley  for  carrier 
qualification  (VT  -  5).  At  this  point  another  split  is  made,  some 
students  proceed  to  Corpus  Christi,  Texas  for  advanced  training  in 
multi-engined  aircraft,  while  others  report  to  Sherman  Field  (VT  -  6) 
for  advanced  instruction  before  beginning  heliocopter  training  at 
Ellyson  Field  (HT  -  8). 
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APPENDIX  D 

Simplified  Example  of  Allocation  Model  for  a  Single  Station 

For  purposes  of  illustration  the  allocation  model  set  forth  in 
Chapter  IV  is  applied  to  a  simplified  air  station  consisting  of: 

1 .  Three  cost  centers: 

a.  Administration  -  Cost  Center  A 

b.  Aircraft  Maintenance  -  Cost  Center  B 

c.  Flight  Training  -  Cost  Center  C 

2.  Four  Variable  Labor  Inputs 

a.  Students  -  x, 

b.  Airmen  -  E-2/E-3  -  x2 

c.  Aviation  mechanics  -  x3 

d.  Flight  instructors  -  x« 

3.  Three  Intermediate  Products 

a.  Students  processed  -  y, 

b.  Repaired  aircraft  -  y2 

c.  Aircraft  requiring  repair  -  y~ 

4.  Two  Final  Outputs 

a.  Student  attrition  -  z, 

b.  Trained  pilots  -  z2 

The  values  of  the  technological  coefficients  (a..,  b..,  and  c.  for 
the  respective  cost  centers)  are  determined  by  the  use  of  material 
balance  equations  relating  inputs  and  outputs  for  each  productive  process 
in  each  cost  center. 
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The  following  flow  diagram  represents  the  station  thus  generated 


Cost  Center 


*2 — H 


For  this  example  only  one  productive  process  considered  in  each 
cost  center,  yielding  the  following  vector  equations  for  each  cost 
center: 
Cost  Center  A: 


an 

a21 

a31 

v 
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Cost  Center  B 


'blll 

1 

|    b21 

j-b31 

!-b41 

i 

1 

1  b5i! 

I    b6ll 

[b7lj 

:z2 

i 

iyi 

I 

i 

*3 

Xr, 


(2) 


Cost  Center  C 


11 


(3) 


^2 
'21  I  y3 

:31 

:41 

where  a,,  b,,  y-i  j  are  the  activity  levels  at  which  the  first  productive 
process  in  each  cost  center  is  to  be  operated.  As  indicated  the  sign 
convention  will  be  used  that  intermediate  goods  produced  are  positive, 
those  consumed  are  negative. 
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Combining  equations   (1),   (2),   and   (3)  yields 


lll 


0 

-b 

0 

b 

a21 

0 

a31 

b 

0 

0 

0 

b 

'11 


21 


'31 


41 


51 


61 


71 


'11 


-c 


21 


'31 


'41 


1 


!     Z 


2 

yl 


1 


B 


H 


I 

y3i 

xi| 

i 

X2  ! 


Note  that  the  partition  of  the  coefficient  matrix  indicated  corresponds 
to: 


.A 


(1) 
(2) 
(3) 


Assuming  constraints  on  inputs  and  outputs  are  set  as 


11 


Z2  -  Kl 


Xl  1  K12 


K2l  -  X2  -  K22 


'31 


'41 


'32 


^4  -  K42 
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The  problem  can  be  written  in  the  standard  form  as 

minimize  (c,   a,,-,  +  c2  a-,)   u,   +  (c2  bg,  +  c,  b 7 .. ) 3 -.   + 

(c2  c31  +  c3  c41)Yl 
subject  to: 


b21        ^      »      K1 


all   al                                            "b31  61  1  ° 

-b41  B]  +  cn  Y]        >  0 

b51  B-,  -  c21  Y]        1  0 

kll  -a21   al  -  k12 

k21  _  a31  a1                                         +b61  ^  +  C3]   y,         1  k22 

k31  -  C41   Yl         -  k32 

k41  1                                                         b71  el  1  k42 


a-i  ,    3i ,    y-|  >   1     0 
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APPENDIX  E 
Simplified  Model  of  Centralized  Planning  for  CNABATRA 

For  illustration  of  the  Centralized  Planning  Model  (CPM)  presented 
in  Chapter  VI,  consider  CNABATRA  to  be  composed  of  two  training  fields, 
NAS-1  and  NAS-2.  Each  of  these  will  be  considered  to  be  structured  much 
the  same  as  the  SSM  example  in  Appendix  B.  That  is,  with  but  three 
cost  centers:  Administration,  Aircraft  Maintenance,  and  Flight  Training 

Let  the  flow  of  products  within  and  through  the  system  be  as  shown 

in  Figure  E-l . 

Cost  Center  A.  -  Administration  at  station  j. 
J 

Cost  Center  B.  -  Flight  training  at  station  j. 

Cost  Center  C.  -  Aircraft  maintenance  at  station  j. 

x,.  -  student  input  to  station  j. 

Xp  •  -  E-2/E-3  input  to  station  j. 

x., .  -  aviation  mechanic  input  to  station  j. 

x4-i  -  flight  instructor  input  to  station  j. 

y-,  .  -  students  processed  by  administration  at  station  j. 

y?.  -  "down"  aircraft  at  station  j. 

y^.  -  repaired  aircraft  at  station  j. 

Z-,  .  -  student  attrition  at  station  j. 

1  J 

Zp2  -  qualified  pilots  (NAS-2  only) 

t-. p  -  partially  trained  pilots  transferred  from  NAS-1  to 
NAS-2  for  advanced  training. 

(Note  that  t,?  is  the  same  product  as  z~  in  the  SSM,  Appendix  D). 
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Assuming,  as  before,  one  technology  per  cost  center,  vector  equations 
for  each  of  the  stations  are  as  follows. 


For  NAS-1: 


Cost  Center  A 


Cost  Center  B 


all" 

a21 

.a31. 

"bn 

b21 

b31 

b4l! 

j 

1 

b5l| 

t 

.  b6lj 

w 


11 


i  x 


n 


[x21 


'12 


!zn 

j  r12 
i 

!Yl1 

ki 

> 

jy3l' 

i         j 

1  X41  : 

Cost  Center  C 


For  NAS-2: 


-c 


11 


21 


'31 


w 


13 


21 


31 


31 


Cost  Center  A, 


-d 


11 


*12l 


'21 


'21 


y12! 


J31 


72 
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Cost  Center  B. 


11 


■21 


'21 


'31 


'22 


^22 

yll 


■41 


-e 


51 


'6Zj 


22i 


32 


x42 


Cost  Center  C?: 


-f 


I    f 


11 
22 


il  f 


31 


'23 


y22l 


y32 
x. 


v32 
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Equations  can  be  grouped  as  before,  yielding: 


!  0 

bll 

0 

0 

0 

0 

jO) 

0 

0 

0 

0 

en 

0 

0 

0 

0 

0 

e21 

0 

a(D* 

0 

b21 

0 

-dll 

0 

0 
0 

all 

"b31 

0 

0 

0 

0 

b41 

-cll 

0 

0 

0 

A<2> 

0 

"b51 

c31 

0 

0 

0 

*wiT 

0 

0 

0 

d21 

"e31 

0 

w12 

0 

0 

0 

0 

G41 

-fll 

0 

0 

0 

0 

"e51 

f21 

w13 

i 

w21 

a21 

0 

0 

0 

0 

0 

a31 

0 

0 

0 

0 

0 

w22 

o 

0 

c31 

0 

0 

0 

w23 

j 

A<3> 

0 

b61 

0 

0 

0 

0 

0 

0 

0 

d31 

0 

0 

0 

0 

0 

0 

0 

f31 

^0 

0 

0 

0 

e61 

0     ' 

y  a(.3)w. 
3-i  J     J 

a21 

0 

0 

0 

0 

0    ~ 

■wir 

a31 
0 

0 
0 

0 
c31 

d31 
0 

0 
0 

0 
f31 

w12 

w13 

.0 

b61 

0 

0 

e61 

0 

w21 
w22 

\A/  _    _ 

z21 

*11 

*21 

^31 
^12 
^22j 

"  i 

xn; 
x2ii 

x   I 
x31  I 

X41  j 

1 

X22| 

x32 

x42 


xl 

x2 

X3| 

1 

X 

'23  j 

The  constraints  imposed  by  Aj  '   are  independent  of  those  imposed  by 
A^]),  and  similarly  for  the  A^  and  fS.3\     j  =  1,2. 

J  J 
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Assuming  the  following  bounds  have  been  specified: 


z22  1  K1 


"   -  x   -  3 
ij     ij     ij 


and, 


i  =  1.....4 
J  -  1.2 


j  =  l1J    ' 


The  problem  is  specified  by 
minimize 


(Cna21  +  C21a31)  wn  +  (C41b6])  w]2  * 


(C^c-J  w1Q  +  (C00d,n)  w01  +  (C.0ecl)  w 


'3r3r  "13 


22U3T   "21 


'42C6T   "22 


(C32f31)  w23 


subject  to: 


allwll 


b21w12 
■b31w12 


-dllw21 


b41w12 


■Cn  nW- 


b51w12  +c31w13 


e21w22 


d21w21  -e31w22 
e41w22 


> 

Kl 

> 

0 

> 

0 

> 

0 

> 

0 

.fnw23l0 


-e51w22     +f21w23i0 


"11  - 

a21 

wll 

"21  - 

a31 

wll 

'31  - 

c31w13 

"41  - 

b61w12 

"22  - 

d31w21 

"32  - 

"42  - 

<  Bn 


<   6 


<   3 


21 


31 


e61w22 


<  Hi 

1622 

f31w23l  632 

ie42 
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and: 


\ 

ia21wll 

lbl 

h 

ia31wll 

+d31w21 

lb2 

h 

< 

c31w12 

f31w23^-  b3 

u 

< 

b61w12 

e61W22              -  b4 

w.  .  =  0,  for  all  i  and  j 


Solution  of  the  L.P.  will  yield  an  optimal  set  of  activity  levels 


w 


w 


w 


11 

* 

12 

• 

w13 
i  * 
w21 

iw*2 
I  * 
w23 


The  optimal  manning  levels  can  be  determined  by 


11 

=         a21wll 

• 

* 

21 

=      a3iwn 

* 

31 

=         c31w31 

* 
41 

=         b61w12 

* 
22 

=         d31w21 

* 
23 

=         f31w21 

• 

• 

24 

=         e61w23 
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APPENDIX  F 
DECOMPOSITION  PRINCIPLE 

The  decomposition  principle  is  a  technique  which  aids  in  the 
optimization  of  a  complex  supersystem  composed  of  many  subsystems  which 
operate  almost  independently  but  have  a  few  constraints  and  an  objective 
function  which  are  held  in  common.  The  supersystem  is  decomposed  into 
subprograms  which  correspond  to  each  of  the  independent  parts,  and  a 
master  program  which  ties  the  subprogram  together. 

The  technique  has  the  disadvantage  that  the  master  and  subprograms 
may  require  solution  several  times. 

The  general  form  of  decomposition  will  be  discussed. 

A.  GENERAL  FORM 

Consider  a  supersystem  made  up  of  two  such  subsystems.  The  linear 
program  for  the  supersystem  may  be  formulated  in  the  following  manner: 

minimize  C,  X  +  C?  Y 

subject  to  A1  X  =  b]  (1) 

A2  Y  =  b2  (2) 

A,  X  +  A2  Y  =  b  (3) 

Where  X  and  Y  are  vectors  of  activity  levels  in  the  first  and  second 
subsystem  respectively.  Equation  (1)  expresses  the  constraints  which 
involve  only  the  first  subsystem,  equation  (2)  is  the  constraints  only  on 
the  second,  and  equation  (3)  expresses  those  constraints  held  in  common. 
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Let  S,  and  Sp  be  the  set  of  extreme  point  solutions  to  A,  X  =  b, 
and  A?  Y  =  b2  respectively.  It  is  assumed  that  the  elements  of  S,  and 
S?  form  closed  bounded  sets.  Under  this  assumption,  any  X  a  feasible 
solution  of  A,  X  =  b,  can  be  represented  as  a  convex  linear  combination 
of  the  elements  of  S, .  Consider  the  elements  of  S,  to  be: 

S-i   =   X-.  ,  Xp »  .  . .  ,  X. 

Then  any  solution  X  can  be  written: 

K 

x     =     y   x.  x. 
3-i  J   J 


where 

K 


}     x.   =   1 ,   x.  >  0 


Similarly  any  feasible  solution  Y  to  Ap  Y  =  bp  can  be  represented 
by  a  convex  combination  of  the  elements  of  Sp.  That  is: 

L 

Y   =   V   y.  Y. 

3-1  J   J 

L 
I      y .   =   1,   y  '  >   0 
j  =  l  J  J 

where, 


S2  =  JYi  •  Y2 YLj 

Thus  for  any  feasible  solutions  X  and  Y,  the  supersystem  program  can  be 
re-written  in  terms  of  the  x.  and  y.,  as  follows: 

T    K         ..    K 

minimize        Z\   (  J  x.  X.)  +  Ci  (  Y  u-  Y.) 

I   j=1  J  J    ^   j=1  J  J 

subject  to  K  L 

At  (I      A  X  )  +  A  (I        y  Y  )    =  b 
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K 


L 


=    l 


=    l 


V  pj 


Conversely  any  a.     and  p.  satisfying  the  above  determine  solutions  X 
and  Y  which  are  feasible  solutions  to  the  supersystem  program. 
Define  the  vectors: 


Ah 

d  x2 


c1  xK 


,     G, 


Ah 

C      Y 
L2  Y2 


cl  L 


and  the  matrices 


[A1   X]f  A]   X2   , 


'  Al   XK] 


'    O  L"p       '   "I    J      "O       •  p       »•••)      r\p       l|J 


Designate  the  columns  of  F,  as  (fi-ii  f-ip> 
Then  the  supersystem  program  becomes: 
minimize 


,  f,K) ,  similarly  for  F«, 


subject  to 


G-i   A   +   Gr,  \i 


F1  A  +  F2  ji  =  b 


K 

I 

3=i 

A.                   =      1 
J 

L 

j-1  J 

\  .,   )i.    ■     0 
J       J  - 
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The  problem  thus  formulated,  with  all  extreme  point  solutions  to  the 
systems  A,  X  =  b,  and  Ap  Y  =  bp  known,  is  termed  the  Extremal  or 
Master  program.  It  has  the  drawback  that  the  set  of  all  extreme  point 
solutions  will  not,  in  general,  be  known  and  will  be  difficult,  (and 
unnecessary)  to  obtain. 

Assume  that  some  basic  feasible  solution  to  the  master  program  is 

known.  A  new  vector  would  be  considered  for  entry  if  its  "Z.  -  C."  is 

J    J 

greater  than  zero.  In  the  Simplex  L.P.  Algorithm,  Z.  -  C.  is  calculated 
from  the  formula  CD  B~  a.  -  C,  where  CD  is  the  transpose  of  the  costs 
of  variables  in  the  basis,  B"   is  the  universe  of  the  present  basis, 

and  a.  is  the  vector  of  coefficients  from  the  A  matrix.  From  Duality, 

T  -1  3 

CB  B    is  equal  to  the  vector  of  dual  variables,  V. 

Consider  that  the  right  hand  side  b,  represents  m  constraints. 

Then  the  dual  of  the  master  program  will  have  m  +  2  variables.  The 

vector  a.  can  be  represented  as: 


:fu 


a.  = 

J 


or 


W 

0 

1 

Thus  for  vectors  from  Subprogram  1 : 


Z.  -  C.  =  Vf,.  -  C.  +  V,, 

J    J      lj    J    m+1 


For  vectors  from  Subprogram  2 


Z.  -  C.  =  Vf,.  -  C.  +  v  ,, 

J    J      2j    j    m+2 


Hadley,  G. ,  Linear  Programming,  p.  229,  Addison-Wesley,  1963. 
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Recall 


lj      1   J 


f0.     =     A9  Y. 
2j      2  j 


In  order  to  choose  a  new  vector  to  enter  the  basis,  the  subsystems 
are  called  upon  to  find  another  extreme  point  solution  to  each  subprogram. 
Each  subprogram  determines  in  effect,  which  of  its  solutions  yields  a 
maximum  "Z.  -  C."  This  is  brought  about  by  the  master  programmer 
specifying  new  "prices"  for  the  X's  and  Y's.  The  new  problems  specified 
to  the  subsystems  are  the  following: 

Subprogram  1 : 

maximize  (VA,  -  C,)  X 

subject  to         A]  X  =  b]  (4) 

X  >  0 

Subprogram  2: 

maximize  (vX,  -  C«)  Y 

subject  to         A2  Y  =  b2  (5) 

Y  >_  0 

Let  the  solutions  to  (4)  and  (5)  be  C,  and  C0       ,.   m3C+.oy,  nvlftnM_,a, 

1      I.     The  master  programmer 

selects: 

min  [C}   +  vm+1,  C2+  vm+2]  >  0  (6) 

and  enters  the  \K+-,  or  p.  .,  which  corresponds  to  (6)  into  the  basis. 
Consider,  for  example,  that  \K+1  will  be  brought  into  the  basis.  The 
new  master  problem  becomes: 


56 


minimize 


V  +  Cl  XK+1  AK+1  +  G2P 


subject  to 


F1A  +  Al  XK+1  AK+1  +  F2P 


K 


j-.i  XJ  +  AK+1 


=  b 
=  1 


A  .,  u.  =  0 
J    J 


Optimality  is  obtained  if  at  some  step: 


rci +  W  <  ° 


and 


tC2  +  W  <  ° 
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